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ABSTRACT 


W49N is a mini-starburst in the Milky Way and thus an ideal laboratory for high- 
mass star formation studies. Due to its large distance (11.kpc), the kinematics 
inside and between the dense molecular clumps in W49N are far from well under¬ 
stood. The SMA observations resolved the continuum emission into two clumps. The 
molecular line observation of S0 2 ( 284 , 24 - 283 , 25 ) suggests that the two clumps have a 
velocity difference of ~7 km s -1 . The eastern clump is very close to two radio sources 
“Gl” and “G2”, and the western clump coincides with a radio source “B”. The HCN 
(3-2) line reveals an extremely energetic outflow, which is among the most energetic 
molecular outflows in the Milky Way. This is the first report of high-velocity molecular 
outflow detection in W49N. The outflow jet might be in precession, which could ac¬ 
count for the distribution, velocity and rotation of water maser spots. Three absorption 
systems are identified in HCO + (3-2) spectra. The absorption features are blueshifted 
with respect to the emission of S0 2 ( 284 , 24 - 28315 ) lines, indicating that a cold layer is 
expanding in front of the warm gas. Further analysis indicates that the expansion is 
decelerated from the geometric expansion centers. 

Subject headings: stars: formation — ISM: kinematics and dynamics — ISM: jets and 


outflows 
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1. Introduction 


High-mass stars play a major role in the evolution of galaxies. When compared with their 


infall rates in 


ligh-mass star formation 


Beuther et al. 


2002; 


Zhang et al. 


regions 


2001 . 


2005 


are orders of magnitude larger ( 

Wu et al. 

2004; 

Oiu et al. 

OO 

o 

o 

(N 

Liu et al. 

2013tj), indicating 


that high-mass star formation is a more energetic process. Additionally, in contrast to 
isolated low-mass star formation, most of high-mass stars form in clusters. However, the 
formation of high-mass stars is still far from clearly understood due to difficulties in observations 


caused by their large distances, clustered environments,jffiort evo 


contamination from the feedback of protostars (Zin neck er & Yorke 


utionary timescales and 


20071) . To study physical and 


chemical environments of high-mass star forming regions, high spatial resolution interferometric 
observations are required. 


Located at a distance of 11.1^ kpc (Zhang et al. 


20131) . W49 is amo ng the most luminou s 


(> 10 7 LA) and massiye (~ 10 6 M 0 ) star forming regions in the Milky Way (ISievers et al. 


Galvan-Madrid et al. 


1991 


20131) . making it an ideal laboratory for high-mass star formation studies. 


As shown in the left panel of Figure 1, W49 contains three massive star forming regions, W49 
north (W49N), W49 south (W49S) and W49 southwest (W49SW), among which W49N is the 
most prominent one containing 40-50 of U C Hn regions which are arranged in a remarkable 2 pc 


diame ter ring (“Welch ring”) (I Welch et al. 


1987 


De Pree, Mehringer & Goss 


1997 


Smith et al. 


2009). The “Welch ring” is shown in the Upper-right panel of Figure 1. 


Several scenarios have been proposed to expl ain the triggering mechanism of th e 


mini-starburst in W49N, including global collapse 

(Welch et al. 

1987; \ 

Villiams et al. 2004 


Galvan-Madrid et al. 

2013b, cloud-cloud collision ( 

Vlufson & Liszt 1977 

; Mivawaki, Havashi, & Hasegawa 

1986; 

Tarter & Welch 

1986; Serabvn, Guesten,& Schulz 1993; 

Bucklev & Ward-Thompson 


1996; 

De Pree, Mehringer & Goss 1997) and expanding shells ( 

Peng et al. 201(]). The global 
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collapse scenario was proposed based on the detection of redshifted 


P-Cygni profile”) in HCO + (1-0) (Welch et al. 


absorption featurcs_(“inverse 


1987) and CS (2-1) (Williams et al. 


20041) lines. 


The Cloud-Cloud collision scenario is mainly based on t 

observed in molecular lines at different radial velocities ( 

te presence of overlapping 

clumps 

Serabvn, Guesten,& Schulz 

1993; 

Bucklev & Ward-Thompson 

1996; 

De Pree. Mehringer & Goss 

1997 

). 

Peng et al. 

201oj) 


identified two expanding shells in the mid-infrared images and confirmed by follow-up molecular 
line observations. They argued that the expanding shells in W49N are suggestive of triggered 


massive star formation in just ~10 5 yrs. Recently, 


Galvan-Madrid et al. 


(2013|) argued that 


the W49N starburst most likely formed from global gravitational contraction with localized 
collapse in a “hubfilament” geometry through large mosaic mapping observations with 
the Submillimeter Array (SMA) in combination with single-dish observations. Previous 
works have indicated that W49N contains substructures and forms a cluster of OB stars. 
However, the kinematics and interactions of different components in this region are still 
far from well understood. In this paper, for the first time, we report the discovery of an 
extremely energetic outflow and localized expansion in the central region of W49N from the 
Submillimeter Array (SMA) observations with an angular resolution of ~ 2".4. 


2. Observations 


The SMA dataset is from the SMALU archive (Ho, Moran, & Lo 


aQ; 


2004). The observations of 


W49N were carried out in 2005 September in its compact configuration. The phase reference 
center was R.A.(J2000) = 19 h 10 m 13.41 s and DEC.(J2000) = 09°06' 14.30". The 345 GHz receivers 


'The Submillimeter Array is a joint project between the Smithsonian Astrophysical Obser¬ 
vatory and the Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the 
Smithsonian Institution and the Academia Sinica. 
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were tuned to 265 GHz for the lower sideband (LSB) and to 275 GHz for the upper sideband 
(USB). Each sideband has 3072 channels and a total bandwidth of ~2 GHz. The frequency 
spacing across the spectral band is 0.812 MHz or ~0.9 km s _1 . In the observations, QSOs 
3C454.3, 1751+096, and planet Uranus were observed for bandpass correction, antenna-based 
gain correction and flux-density calibration, respectively. 


MIRIAD was employed for calibration and imaging (Saultet al 


1995 ). The 1.1 mm 


continuum data were acquired by averaging all the line-free channels over both the upper and 
lower spectral bands. MIRIAD task ’’selfcal” was employed to perform self-calibration on the 
continuum data. The gain solutions from the self-calibration were applied to the line data. 


The primary beam of the SMA at 270 GHz is ~40". The uv distances range from 10 

to 62 kA. The synthesized beam size and 1 cr rms of the continuum emission are 2.77" x 2.02" 
(PA=88.9°) and ~20 mJy beam -1 , respectively. We smoothed the spectral lines to a spectral 
resolution of 1 km s -1 and the corresponding lcr rms for lines is ~0.1 Jy beam -1 per channel. 


3. Results 

3.1. 1.1 mm continuum emission 

The 1.1 mm continuum image is shown in red contours and color image in the lower-right 
panel of Figure 1. Two clumps are revealed in the 1.1 mm continuum emission. The 1.1 mm 
continuum emission can be well fitted with two 2-dimensional Gaussian sources as shown in blue 
contours. The peak positions of the eastern and western clumps are (-0".35+0".09, -1". 13+0".03) 
and (-4".06+0". 13, -1".09+0".06) offset from the phase center, respectively. The total integrated 
flux and deconvolved size of the eastern clump are 4.40 Jy and 2".04 x 0".71 (P.A.=-74.7°), 
respectively. Assuming a distance of 11.1 kpc, the effective radius of the eastern clump which 
is defined as ^^ab, where a and b are the full width at half-maximum (FWHM) deconvolved 
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sizes of major and minor axes, is ~0.065 pc. The western clump is not resolved by the SMA 
and appears to be a point source with a total integrate d flux of 2,02 Jy. In the upper pan el of 


Figure 2, we present the 3.6 cm continuum emission (De Pree, Mehringer & Goss 


1997) as a 


color image and also label the previously discovered bright radio continuum sources. We find the 
eastern continuum emission peak is very close to two radio sources “Gl” and “G2”. The western 
continuum emission peak coincides with a radio source “B”. 


3.2. Line emission 

Dozens of molecular lines were detected in both sidebands and were identified with the 
Splatalogue database for astronomical spectroscopy. Here we report the results from S0 2 
(28 4 ,24-28 3 , 2 5), HCN (3-2), and HCO + (3-2) lines. 


3.2.1. S0 2 (2§424-28 3 25 ) emission 

The integrated intensity from -10 to 30 km s -1 of S0 2 ( 284 , 24 - 283 , 25 ) is shown in gray contours 
in the upper panel and black contours in the lower panel of Figure 2, respectively. The emission 
of S0 2 (28 4 , 24 - 283 , 25 ) coincides with the 1.1 mm continuum emission very well. From its intensity 
weighted velocity map (moment 1) in the lower panel of Figure 2, one can see that there is a large 
velocity gradient from east to west. This large velocity gradient is more clearly seen from the 
position-velocity diagram of S0 2 ( 284 , 24 - 283 , 25 ) (Figure 3), which is cut along the black dashed 
line in the lower panel of Figure 2. The peak velocities of S0 2 ( 284 , 24 - 283 , 25 ) are estimated from 
Gaussian fits at the individual positions and are marked with blue crosses in the lower panel of 
Figure 3. The peak velocities of S0 2 ( 284 , 24 - 283 , 25 ) at the continuum emission peaks of the eastern 
and western clumps are 5.8+0.1 and 12.7+0.1 km s -1 , respectively. The velocity gradient between 
these two clumps is ~34.7 km s ^ 1 pc -1 . 
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3.2.2. Extremely high velocity emission in HCN (3-2) lines 


As shown in the position-velocity diagram (upper panel of Figure 3) and averaged spectrum 
(upper panel of Figure 4), the HCN (3-2) line emission shows extremely high velocity emission. 
The terminal velocities of the blueshifted and redshifted emission are as high as ~-70 and ~90 
km s _1 , respectively. The redshifted high-velocity emission is contaminated by HCN v2=l (3-2) 
line, located at around 50 km s -1 as marked by the dashed dark brown line in the Position-Velocity 
diagram. There is also an emission peak in the averaged spectrum of HCN (3-2) (upper panel 
of Figure 4) at around 50 km s _I due to the contamination of HCN v2=l (3-2) line. To confirm 
the contamination of HCN v2=l (3-2), We integrated emission from 49 to 51 km s _1 and 
present the integrated intensity in Figure 5 as color image and white contours. The spatial 
distribution of HCN v2=l (3-2) emission, even though being contaminated by HCN (3-2) 
high velocity emission, is roughly consistent with 1.1 mm continuum but significantly 
different from outflow emission as shown in Figure 4. The integrated intensity map in Figure 
5 also rules out the possibility of contamination emission from foreground clouds, which 
might show a very extended structure if exists. 


Roberts et al. 


(2011) observed the J=3-2 and 4-3 transitions of HCN with the JCMT 


toward W49N. However, due to poor angular resolution and sensitivity, they did not detect 
such high velocity wings. Their spectral can be well fitted with two gaussian components, 
indicating the existence of substructures. The HCN outflow discovered in this work should 
be the first report of high velocity outflow detection in W49N. We divide the outflow emission 
into “Low” and “High” components. The velocity intervals for each component are shown in the 
second column of Table 1. The emission between 40 and 60 km s _I is not included in further 
analysis due to the contamination of HCN v2=l (3-2) line. The integrated intensity maps of 
the blueshifted and redshifted “High” outflow emission are shown as blue and red contours, 
respectively, in the lower panel of Figure 4. The blueshifted and redshifted outflow emission 
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overlaps and only displaces by 0".7 (AO) in projection, indicating the outflow axis is almost 
along the line of sight. The projection deconvolved angular size of minor axes of outflow lobes 
is about 2" ( 0 ). The upper limit of the outflow inclination angle with respect to the line of 
sight is y = y = ^ x - 20°. The characteristic outflow velocity (V c har) is defined as 
Vchar = V flow - V sys , where \ sys is the systemic velocity (5.8 km s' 1 ) of the eastern continuum 
clump and V f Iow is the intensity weighted velocity of high velocity emission corrected with the 
projection effect by assuming an outflow inclination angle of 20°. The characteristic outflow 
velocity of each outflow component is presented in the third column of Table 1. The V c har of the 
high velocity redshifted outflow (“High Red”) component is as high as 71.5 km s _1 , indicating 
that the outflow is extremely energetic. The effective emission radius of each outflow component 
shown in the fourth column of Table 1 is derived from its deconvolved size and also corrected with 
the projection effect. The effective radii of the emission areas of the “Low” and “High” outflow 
components are 0.1 and 0.2 pc, respectively. The outflow properties will be discussed in detail in 
section 4.1. 


3.2.3. Absorption in HCO + (3-2) lines 

Figure 6 presents spectra of HCO + (3-2) and S0 2 ( 284 , 24 - 283 , 25 ) at six positions. Three 
absorption systems in HCO + (3-2) lines are identified. The mean velocities of these three 
absorption systems are “-7.5 km s _1 ”, “2.0 km s _1 ”, and “12.5 km s _1 ”, respectively. These three 
absorption systems can also be clearly seen from the Position-Velocity diagrams in Figure 3 as 
in dashed contours. We made the integrated intensity maps in the velocity intervals of [-10,-4], 
[-2,6] and [10,15] km s _1 for the “-7.5 km s _1 ”, “2.0 km s _1 ”, and “12.5 km s -1 ” absorption 
systems, respectively, and display them as dashed contours in the upper panel of Figure 2. The 
“-7.5 km s^ 1 ” absorption is associated with radio continuum sources “H” and “G4”. The “-2.0 
km s _1 ” absorption perfectly corresponds to the 3.6 cm continuum emission of “G”, which 
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includes sources “Gl” to “G5”. The “12.5 km s _1 ” absorption is associated with radio source “B”. 
Additionally, from the spectra in Figure 6 and Position-Velocity diagram in Figure 3, one can see 

that HCO + (3-2) absorption dips are always at the blueshifted side of the S0 2 (28 4 i 24 -28 3i25 ) 
emission except in the overlapped regions (e.g. (-1".5,-1".5) and (-2".5,-1".5) off positions) 
of the two clouds. The HCO + (3-2) spectra in the right column of Figure 5 have a typical 
“P-Cygni” profile with emission peaks redshifted with respect to absorption dips. The peak 
velocities of HCO + (3-2) absorption dips estimated from gaussian fit at each position are also 
marked with red crosses in the lower panel of Figure 3. One can see that the absorption dips 
of HCO + (3-2) are mostly blueshifted to the emission peaks of S0 2 (28 4j24 -2 8 3 , 2 5 ). Since S0 2 
(28 4 24 -28 3 , 25 ) (E„ ~ 416K) has much higher upper energy level than that of HCO + (3-2) 

(E„ ~ 26K), this velocity placement indicates that HCO + (3-2) absorption traces an localized 
expanding cold layer in front of the warm gas traced by S0 2 (28 4 24 -28 3 25 ) emission. More 
discussion about the absorption in HCO + (3-2) lines will be presented in section 4.2. 


4. Discussion 


4.1. Extremely energetic molecular outflow 

4.1.1. The outflow properties 


As presented in section 3.2.2, the HCN (3-2) lines show extremely high velocity emission, 


indicating energetic outflows. Following 
given by: 


Liu et al. 


( 201 lbh . the total mass of each outflow lobe is 


() pEn/Trot f ' 

M = 1.04 x 10 -*D- ot ~ , -- S v dv (1) 

XV S p- J 1 - e T 

where M, D, S y , x, and r are the outflow gas mass in M G , source distance in kpc, line flux density 
in Jy, fractional abundance with respect to H 2 , and optical depth, respectively. Q rot , E„, T rof , v, 
S, p are the rotational partition function, the upper level energy in K, the rotation temperature, 
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the rest frequency, line strength and the permanent dipole moment, respectively. Owing to the 
lack of a direct estimation of the fractional abundance of HCN in the outflow region of W49N, 


2011b) 


we adopt the same fractional abundance of 5.2 x 10 9 as that in G9.62+0.19 F (iLiu et al. 
in calculating the outflow parameters. Assuming that the high-velocity emission is optically thin 


and T, y ,,=30 K (Liu et al 


201 lb), the outflow masses were calculated and presented in the sixth 


column of Table 1. The total outflow mass is ~40 M 0 . The momentum can be calculated by 
P = M x V char , and the energy by E = \M x V^ har - The momentum and energy of the outflows are 
shown in the seventh and eighth columns, respectively. The total momentum and energy are 1700 
M g km s -1 and 8.7 x 10 47 erg, respectively. The dynamical timescale t is estimated as 2R e ff/V char 
and presented in the fifth column of Table 1. The dynamic timescale is about 0.5 - 1 x J_0 4 year. 


In comparison with surveys of molecular outflows in massive star forming regions (Zh ang et al. 


2001 . 


2005), the outflow mass, momentum and energy in W49N are in the upper bound of the 


outflow parameters of massive star formation. 

The mechanical luminosity L„„ and the mass-loss rate M are calculated as L,„=E/t 
and M = P/(tV w ). Here we assume that the outflow is powered by winds driven by 


accretion disks and the wind velocity V H , is 500 km 1 (Zhang et al 


2005). The mechanical 


luminosity and the mass-loss rates are listed in the last two columns of Table 1. The total 
mechanical luminosity and the mass-loss rate of the outflows are estimated to be 1.3xl0 3 L 0 and 
5.6 x 10“ 4 Mo-yr^ 1 , respectively. Assuming that only one third of accreted mass is ejected in winds 


(Shu, Adam & Lizano 


1987), we find a lower limit to the accretion rate to be 1.7 x 10 3 M 0 -yr 1 . 


Smith et al. I (2009) identified a mid-IR source (G:IRS1) and assigned it to the driving source 


of the outflow. From SED modeling, they inferred that the embedded protostar has a stellar mass 
of ~45 M 0 , luminosity ~ 3 x 1O 5 L 0 and an equivalent envelope accretion rate ~ 10 3 M 0 yr 1 . The 
mechanical luminosity of O-type stars with a mass loss rate of M and stellar wind velocity v w can 
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be estimated from (INagv et al. 1120121) : 


M 


L w = 1.3 x 10 ib (- £ 

10 6 M q yr 1 2 x 10 3 km s 


t)(- 




-ferg s" 1 


( 2 ) 


The mechanical luminosity of an O type star with a mass loss rate of 10 6 M G yr 1 and stellar wind 


velocity of 2000 km s 1 (Smith 


2006 ) is ~340 L 0 . Assuming an efficiency of 10%, its mechanical 


heating rate (34 L G ) is much smaller than the mechanical luminosity (1.3xl0 3 L 0 ) of outflows, 


indicating that the outflow could not be driven by the stellar winds o 
G:IRS1). The outflow might be driven by accretion as suggested by 


a single O 


t ype star (e.g. 


Smith et al. 


enve 


(2009). The 


ope accretion rate is in their model fits to the IR SED is 1 x 10 3 M G yr 1 (iSmith et al. 
2009), which is consistent with the accretion rate (1.7 x 10 -3 M G yr 1 ) inferred from the HCN 


outflow. 

There are large uncertainties in the estimation of the outflow parameters. The optically thin 
assumption may not be true for low-velocity outflows. Assuming an average optical depth of 
1.5, the corresponding outflow parameters (M, P, E, M and L m ) will become two times larger. 
As mentioned before, we ignored the emission at velocities between 40 to 60 km s _1 due to the 
contamination of HCN v2=l (3-2) line, which definitely leads to underestimation of the outflow 
parameters for redshifted high-velocity emission. The outflow parameters are less sensitive to 
the temperature. If we take a rotation temperature of 60 K (two times larger than the value we 
used), the outflow parameters only become 30% smaller. The largest uncertainties come from 


the fractional abundance of HCN and inclination angle. 


Roberts et al. 


(2011) derived an Hi 


volume density n ff , of 2xl0 6 cm 3 and an HCN column density of (2-7)xl0 15 cm -2 in the 
“Core” region of W49N. The angular size of HC 15 N (4-3) emission is L ~ 23". We can 
estimate an average H 2 column density: N Hl = Lx N Hl = 7.6 x 10 24 cm 2 . Therefore, the 
fraction abundance of HCN at “Core” region is (2.6 - 9.2) x 10 _1 °, which is about one 
order of magnitude smaller than the value (5.2 x 10 -9 ) we used, indicating that the outflow 
parameters could be significantly underestimated. However, HCN appears greatly enhanced 
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in some outflow regions. In NGC 1333 IRAS 2A outflow region, thejractional abundance of 


HCN is about twice larger than that in the core region (Jorgensen et al. 


2004 . While in LI 157 


outflow region, the HCN abundance is even enhanced by a factor of 100. Thus, more accurate 
determination of the HCN abundance in outflow regions is critically important in estimation of 


outflow parameters. 


Su et al. 


(2007) estimated an HCN abundance of 1 - 2 x 10 8 in the massive 


outflow lobes of IRAS 20126+4104. If we assume that HCN fractional abundance in W49N is 
lxlO -8 as in IRAS 20126+4104 outflows, the outflow parameters will become about two times 
smaller. The inclination angle directly affects the estimation of R e ff and V r/M) -. The outflow 
parameters of M, P, and E do not change too much (< 10%) when the inclination angle changes 
from 10° to 30°. While the corresponding M and L m can change by a factor of 2. In general, the 
respective uncertainty of outflow parameters due to the effects discussed above could be at least 


one order of magnitude. The outflow in W49N was not detected in CO (2-1) line emission 

13). 


with an angular resolution of 3" (Galvan-Madrid et al. 


2013!). The reason may be that 


the outflow is very compact and warm, which cannot be revealed with low excitation lines. 
Further higher angular resolution and higher sensitivity observations of multiple higher 
excitation transitions of HCN and other molecular lines (e.g. CO, SiO) may help to better 
constrain the outflow properties. 


forming regions 

(Su et al. 

2007 

Liu, Ho & Zhang 

2010; 

Liu et al. 

201 la.1 

+ Zhu, Zhao & Wright 

2011; 

Liu et al. 

2013a 

). Comparing with those outflows, W49N has much larger outflow velocity. 


Among the previous detected HCN outflows, the one in G9.62+0.19 F might be the most energetic 


one (Liu et al. 


2011b). The outflow in W49N has similar mass and momentum to those of 


G9.62+0.19 F, but the energy and mechanical luminosity are about 10 times larger. 


Wu et al. 


(2004) cataloged 391 outflows, most of which were detected by single dishes. When compared 


with those outflows, the one in W49N is also among the most luminous and energetic outflows. 


The maximum energy and mechanical luminosity of the outflows in 


Wu et al. 


(2004) are around 
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10 48 erg and 10 3 L G , respectively, which are comparable to that in W49N. Thus, we argue that 
the outflow in W49N is among the most energetic outflows in the Milky Way. However, we 
should caution that the HCN outflow is barely resolved by the SMA and may contain 
multiple outflows as seen in other high-mass cluster forming regions (e.g. NGC7538S; 


Naranjo-Romero et al. 


(2012)). Therefore, the outflow parameters estimated here may be 


the sum of multiple outflows. 


4.1.2. Does the outflow excite the water masers? 


W49N harbors the most luminous water maser activities in the Milky Way ( Walk er et al. 


1982 

). The water masers in W49I 

4 are mainly located near the 

dn regions “Gl” anc 

“G2” 

(Moran et al. 

1973; 

Walker et al. 

1982; 

Gwinn. Moran & Reid 

1992; 

De Pree et al. 

2000; 


Mcgrath, Gross, & De Pree 


2004). 


Mcgrath, Gross. & De Pree 


(2004) identified 316 water maser 


components with a velocity spanning from -352.1 to 375.5 km s _1 . In the lower panel of Figure 4, 
we overlay the water maser spots on the integrated intensity contours of HCN (3-2) high-velocity 
outflow emission. The water maser spots near “Gl” and “G2” are distributed in three distinct 
clusters. Those maser spots with blueshifted and redshifted high-velocity features are located 
to the west and to the east of the mid-IR source “G:IRS1”, respectively. The tightest cluster of 
maser spots is close to “G:IRS1” and has a linear arrangement elongated along a north-south 


axis. 


Honma et al. 


(2004) detected a strong maser outburst on an arc-like structure in the central 
cluster, which may be due to a shock phenomenon powered by a forming star or a star cluster. 


Recently, from proper motion measurements, 


Zhang et al. 


(2013) found that those water maser 


spots appear to be in expansion and rotation with expanding and rotating velocities of 10.99+1.33 
and 5.28+0.80 km s _1 . 


The water maser emission in W49N has been a puzzle to astronomers for decades. 


Tarter & Welch 


(1986) suggested a cloud-cloud collision model for water maser excitation in 
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W49N. An alternative promising ex planation is that water masers form at the 


cavity walls (or jet-driven cocoons) (Mac Low et al. 


1994; 


Smith et al. 


edge of the outflow 


20090 . Such model for 


water maser excitation requires a shock velocity on the cavity walls of ~20-200 km s _1 . Since 
the outflow velocity of HCN (3-2) is up to ~70 km s _1 , water masers should be easily excited 
at the edge of the outflow cavity walls. Additionally, the low inclination angle of the outflow 
might account for the observed high velocity features of water masers. As mentioned before, 
the blueshifted and redshifted high-velocity emission of HCN (3-2) overlaps and only displaces 
by ~ 0".7 in projection, which is not resolved by the SMA. However, the orientation of the 
HCN outflow is roughly consistent with the water maser outflow. The western cluster of water 
masers is located at the peak of the blueshifted high-velocity HCN emission. While the eastern 


cluster of water maser e mission is mainly at t 


outflow was observed by 


Scoville et al. 


le red lobe of HCN outflow. A ~25 km s 1 CO 


( 1986) . where the blueshifted CO emission is towards 


the east and the redshifted CO emission is to the west. The orientation for the low-velocity CO 
outflow is opposite to the high-velocity HCN outflow, indicating that the outflow jet might be in 
precession. The jet precession naturally explains why the high- and low velocity features coexist 
in the centeral cluster of water masers, while at the western and eastern clusters of waters masers, 


only high-velocity features are observed (IMcgrath. Gross . & De Pree 


2004j) . The expansion of 


the water maser spots might be caused by outflows. The jet precession can naturally explain 
why the maser spots are also in rotation in spite of expansion. However, we should point out 
that the jet precession scenario needs to be tested by higher angular resolution observations 
because: 1. The outflow in W49N is not obvious in new CO (2-1) line observations with the 


SMA 


Galvan-Madrid et al. 


(2013); 2. The HCN outflow is barely resolved by the SMA and 


may contain multiple outflows. 
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4.2. Decelerated expansion? 

As discussed in section 3.2.3, the absorption dips of HCO + (3-2) are all blueshifted with 
respect to the emission peaks of S0 2 ( 284 , 24 - 283 , 25 ), indicating that the cold layers may be in 
expansion. Since there is no S0 2 (284 24-283 25 ) emission corresponding to the “-7.5 km s -1 ” 
absorption system, here we only focus on the “2.0 km s _1 ” and “12.5 km s _1 ” absorption systems. 
The systemic velocities of the warm layers can be well determined from S0 2 ( 284 , 24 - 283 , 25 ) 
emission lines. We define the line-of-sight (los) expansion velocity as V^=V em ,-V a ^, where 

\ emi is the velocity of emission peak of S0 2 ( 284 , 24 - 283 , 25 ) and \ abs is the velocity of absorption 
dip of HCO + (3-2). \ emi and \ abs were derived from gaussian fits and marked with blue and 
red crosses in the Position-Velocity diagram in the lower panel of Figure 3, respectively. 

We derive at different positions along the black dashed line in the lower panel of Figure 2 
and plot versus the RA offsets in the upper panel of Figure 7. We ignored the positions 

where the “2.0 km s 1 ” and “12.5 km s 1 ” absorption systems overlaps such as off positions 
(-1".5,-1".5) and (-2".5,-l".5). 

There are two peaks in the “V^ vs. RA” plot in the upper panel of Figure 7: the 

eastern one around 2".5 and the western one around - 6 ". This “V^ vs. RA” curve can not be 
simply fitted with a spherically expanding model with a constant expansion velocity, in which 
the Position-Velocity curve would be elliptical like. Instead, we propose a toy model, which is 
shown in the lower panel of Figure 7, to explain the absorption of HCO + (3-2). Relative to the 
warm gas layer which is distributed in the central plane, the two parallel elliptical cold layers are 
moving outward. The positions of the cold layer can be described with two parameters: height 
from the central plane (“h”) and distance (“r”) to the geometric expansion center “O” whose RA 
offset is “X0”. The maximum of “h” and “r” are “b” and “a” respectively. Assuming this simple 
geometric model, we find the “V^, vs. RA” curve can be well fitted with Y exp oc r 1 , indicating 
that the expansion decreases with the distance to the geometric expansion center. \ exp is the real 
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y!os 

expanding velocity: \ exp = where 6 is the angle of velocity vector with respect to the 
line of sight. The maximum expansion velocity Y max corresponds to the minimum distance “b”. 
The best fits with smallest;^ 2 are shown in red and blue solid lines in the upper panel of Figure 
7. The parameters [b, a, XO, Y max ] of the best fits to the western and eastern expansion systems 
are [1.18", 2.84", -6.24", 3.37 km s" 1 ] and [1.43", 4.57", 2.52", 5.27 km s -1 ], respectively. 


The dynamical timescale can be estimated as t exp = tJ—. The dynamical timescale for the 
western and eastern expansion systems are 1.8xl0 4 and 1.4xl0 4 yr, respectively. The actual 
masses (M) of the expanding cloud can be roughly estimated from their virial masses M VI > 


(Kauffmann, Pillai, & Goldsmith 


,20131: 


M 


10 3 M ffl 


M vij I c 


l-2<7^rr> 2 <— )l«c 


km s 


pc 


(3) 


where a cr is the so-called critical virial parameter and is found to be 2±J_for a wide range of cloud 


2013) 


shapes and density gradients. Here we take a cr of 2 (Kauffmann. Pillai. & Goldsmith 

The radius R e ff here is taken as a effective radius of 'Jab. The R e ff of the western and eastern 
expanding clouds are ~0.10 and —0.14 pc,respectively. Velocity dispersion cr v is measured from 
the linewidth of the absorption lines of HCO + (3-2). The cr v is found to be 1.5 and 1.8 km s _1 
for the western and eastern expanding clouds, respectively. Thus, the masses of the western and 
eastern expanding clouds are ~270 and ~ 130 M 0 , respectively. For “subcritical” clouds, which 
are unbound and expanding, the virial parameter should be larger than the critical virial parameter. 
Thus, the actual mass estimated from equation (3) should be taken as an upper limit. We define 
an effective expanding velocity as Y eff =-^- x V max . V e ff of the western and eastern expanding 
clouds are ~2.2 and ~2.9 km -1 , respectively. Then we can estimate the momentum and kinetic 
energy of the expanding clouds as P=MxV e // and E kin = f MxV 2 ^, respectively. The momentum 
of the western and eastern expanding clouds are ~290 and 790 M Q km s -1 , respectively. The 
kinetic energy of the western and eastern expanding clouds are ~ 6.3 x 10 45 and ~ 2.3 x 10 46 erg, 
respectively. The corresponding mechanical luminosity L m is calculated as L, n =E kin /t exp . The L,„ 
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of the western and eastern expanding clouds are 2.9 and 13.7 L 0 , respectively. 


Smith et al. ( 2009!) suggests that the cavity encircled by radio source “Gl” to “G5” is caused 
by the outflow. Assuming an efficiency of 1-10%, the mechanical heating rate of the outflow is 
about 13 to 130 L 0 , which is comparable or much larger than the mechanical luminosity of the 
eastern expanding cloud, indicating that the outflow has sufficient energy to blow the cavity. The 
extent and timescale of the outflow are also comparable to those of the eastern expanding cloud 
which is associated with the cavity. All indicate that the cavity is formed due to the outflow. 
However, this scenario is doubtful because that the outflow appears to have a very small 
inclination angle and the outflow lobes are very compact (0.1-0.2 pc in radius). The red 
outflow lobe is about 2.7" (-0.15 pc) away from the expanding center of the cavity. This 
displacement is comparable to the radii of outflow lobes and could be resolved by present 
SMA beam size (-2.4"). On the other hand, we also notice that the stellar wind of an O-type star 
also has sufficient energy to create such cavity. The mechanical luminosity of an O type star with 
a mass loss rate of 1O~ 6 M 0 yr _1 and stellar wind velocity of 2000 km s^ 1 is -340 L 0 . Assuming 
an efficiency of 10%, its mechanical heating rate is larger than the mechanical luminosity (13.7 
L 0 ) of the eastern expanding cloud. Therefore, we suggest that the expansion of the cavity is 
more likely due to the mechanical heating produced by the stellar winds rather than outflow. The 
western expanding cloud is associated with a radio source “B”, which contains one 06.5 and one 


05.5 ZAMS type stars (De Pree et al. 


20001) . The HCO + lines show typical “P-Cygni” profile 


toward radio source “B”, which indicates the expansion of an Hii region. 


Two expanding shells with an average radius of -3.3 pc have been identified the 
mid-infrared images and the 13 CO J=2-l and C ls O J=2-l data from the IRAM 30-m 


observations (Peng et al 


2010). The average mass of the expanding shells is 2 x 10 4 M G 


with a kinetic energy of - 10 49 erg, which are 2 or 3 orders of magnitude larger than the 
expanding systems reported in this paper. The expansion traced by HCO + (3-2) line is more 
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localized when compared with those infrared expanding shells. The connection between 
HCO + expansion systems and infrared expanding shells are not clear. 


5. Summary 

We have studied W49N with the SMA data in the 1.1 mm continuum and molecular line 
emission. The main results of this study are as follows: 

1. Two clumps are revealed in the 1.1 mm continuum emission. The eastern one is very close 
to radio source “Gl” and “G2”, and the western clump coincides with the radio source “B”. 

2. HCN (3-2) lines show extremely high-velocity emission, indicating energetic outflow 
motions. This is the first report of high-velocity outflow detection in W49N. The total 
outflow mass is 40 M Q . The total momentum and energy are 1700 M Q km s _1 and 8.7 x 10 47 
erg, respectively. The dynamic timescale of the outflow is about 0.5 - 1 x 10 4 year. 

The total mechanical luminosity and the mass-loss rate are estimated to be 1.3xl0 3 L 0 and 
5.6x 10 -4 M 0 yr _1 , respectively. The corresponding accretion rate is at leat 1.7x 10 3 M 0 -yr _1 . 
The orientation for the high-velocity HCN outflow is opposite to the low-velocity CO outflow, 
indicating that the outflow jet might be in precession. The jet precession can naturally explain the 
distribution, velocity and rotation of water maser spots. 

3. Three absorption systems are identified in HCO + (3-2) spectra, which are at “-7.5 km s _1 ”, 
“2.0 km s- 1 ”, and “12.5 km s _1 ”, respectively. The “-7.5 km s _1 ” absorption is associated with 
radio continuum sources “H” and “G4”. The “2.0 km s _1 ” absorption perfectly corresponds to the 
3.6 cm continuum emission of “G”, which includes sources “Gl” to “G5”. The “12.5 km s _1 ” 
absorption is associated with radio source “B”. The absorption of HCO + (3-2) is blueshifted with 
respect to the emission of S0 2 (28 4 ,24-28 3 15 ) lines, indicating a cold layer is expanding in front of 
the warm gas. Analysis toward the “2.0 km s _1 ”, and “12.5 km s _1 ” absorption system indicates 
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that the expansion velocity linearly decreases with the distance to the geometric expansion center. 
The expansion traced by HCO + (3-2) absorption features might be caused by mechanical heating 
of stellar winds from associated O type stars. 
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19 h 10 m 1 4 s .0 13 s . 6 13 s .2 12 s .8 

RA (J2000) 


Fig. 1.— Left: The 3.6 cm continuum emission (Dc Pree. Mehringer & Goss 


1997) is shown in 


contours overlayed on the 3 color composite image from Spitzer/IRAC (3.6 //m in blue, 4.5 /mi in 
green and 8.0 /urn in red). The contours start at the 3cr level (1.0 mJy beam” 1 ). Subsequent positive 
contours are at 2, 4, 8, 16, 32 times the 3cr level. Upper-right: The ring of UCHu regions (“Welch 
ring”). The color image and contours represent the 3.6 cm continuum emission. The names of 
strongest UCHii regions are labeled. Lower-right: The 1.1 mm continuum emission observed by 
the SMA is shown in color image and red contours. The blue contours represent the 2-dimensional 
gaussian fits. The contours are (5,20,50,100,150)x0.02 Jy beam” 1 (1 cr). The SMA beam is shown 
in filled black ellipse. 
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Fig. 2.— Upper: The 3.6 cm continuum emission (Dc Prcc. M eh ringer & Goss 


1997) is shown 


in color image. The names of bright cm continuum sources are labeled from “A” to “H”. The 
integrated intensity maps of “-7.5 km s -1 ”, “2.0 km s -1 ”, and “12.5 km s _1 ” absorption of HCO + 
(3-2) are shown in green, blue and red dashed contours, respectively. The 1.1 mm continuum emis¬ 
sion and integrated intensity map of S0 2 (284 24-283 25 ) are shown in solid black and gray contours, 
respectively. All the contours are started from 30% in step of 20% of the peak values. The peak val¬ 
ues for the green, blue, red, gray and black contours are -12.8, -19.8, -8.7, 106.2 Jy beam -1 km s -1 
and 3.2 Jy beam -1 , respectively. Lower: The intensity weighted velocity (Moment 1) map of S0 2 
(28 4 i 24 -283 ;25 ) is shown in color image. The integrated intensity map of S0 2 (284,24-283,25) is shown 
in black contours. The contours are [-3, 3, 5, 10, 20, 40, 80, 160, 240, 320]x0.3 Jy beam km s -1 (3 
cr). The 1.1 mm continuum emission is shown in blue contours. The SMA beam is shown in filled 
black ellipse. 














-26- 



l 20 -10 0 10 20 


Velocity (km/s) 

Fig. 3.— Position-Velocity diagrams of HCN (4-3), S0 2 (28 4 ,2 4 -28 3 j 2 5 ) and HCO + (3-2) are shown 
in black, blue and red contours, respectively. The Position-Velocity diagrams are made by cutting 
along the black dashed line in Figure 2. The R.A. offsets are relative to 19 h 10 m 13.41 s . All the 
contours in the upper panel are started from 1 Jy beam -1 (10 cr) in step of 1 Jy beam -1 . All the 
contours in the lower panel are started from 0.3 (3 cr), 1 Jy beam -1 in step of 1 Jy beam -1 . The 
dashed straight lines in the upper panel represent the contaminated molecular lines (S0 2 (13 3 j h- 
13 2 ,i2)> OCS (22-21) & HCN v2=l (3-2) ). These dashed straight lines are oblique due to 
velocity gradient. The red and blue crosses in the lower panel mark the peak velocities of S0 2 
(28 4 ,24-28 3 ,2 5 ) emission lines and HCO + (3-2) absorption lines at each position from Gaussian fits. 
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Fig. 4.— Upper: Spectrum of HCN (3-2) line, which is averaged over 4" x 4" area centered 
at position (-0".5,-l".5) offset from the phase reference center. Lower: The 7 mm continuum 


emission (De Pree et al. 


20Q0) is shown in color image. The integrated intensity maps of high 


velocity blueshifted and redshifted emission of HCN (3-2) are shown in blue and red contours. 
The red crosses and blue boxes mark the positions of water masers with positive and negative LSR 
velocities, respectively. All the contours are started from 20% in step of 20% of the peak values. 
The peak values for the blue and red contours are 30, and 49.4 Jy beam -1 km s -1 , respectively. The 
mid-IR source “G:IRS1” is marked with yellow filled star. 
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Fig. 5.— The integrated intensity of HCN v2=l (3-2) line is shown as white contours and 
color image. The 1.1 mm continuum emission is shown in black contours. All the contours are 
started from 20% in step of 20% of the peak values. The peak value of HCN v2=l (3-2) is 13.9 
Jy beam" 1 km s _1 . 










-29- 



Fig. 6 .— Spectra of HCO + (3-2) and S0 2 ( 284 , 24 - 283 , 25 ) at six positions (labeled at the upper- 
right comer of each panel in unit of arcsec, and take along the line shown in Figure2-bottom. 
) are shown in red and blue, respectively. The velocities of the absorption dips of HCO + (3-2) 
and emission peaks of HCO + (3-2) and S0 2 ( 284 , 24 - 283 , 25 ) are marked with red and blue dashed 
vertical lines, respectively. 
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Fig. 7.— Lower: A cartoon of one possible geometric model for causing the absorption in HCO + 
(3-2) lines in W49N. There are two expanding systems. The cold gas traced by HCO + ( 3 - 2 ) 
absorption is distributed in two parallel expanding elliptical layers. The warm gas traced by 
S 0 2 (284,24-283,25) is distributed along the middle plane. The color change of the horizon axis 
reflects the velocity gradient in S 0 2 ( 284 , 24 - 283 , 25 ) emission. The maximum height (“h”) of the 
cold layer is “b”, and the maximum distance (“r”) of the cold gas to the geometric center (“O”) of 
the warm layer is “a”. The RA offset of the geometric center is “X0”. The maximum expanding 
velocity is “V max ”. Upper: The observed line-of-sight (los) expanding velocities at individual 
positions are plotted as crosses with errorbars. The line-of-sight (los) expansion velocity is 
defined as V^=V em! -V ai)i , where V emi is the velocity of emission peak of S0 2 ( 284 , 24 - 283 , 25 ) and 
Y a bs is the velocity of absorption dip of HCO + (3-2). The two 1.1 mm continuum emission peaks 
are shown in stars. The best fits to the position-velocity curves are shown in red and blue solid 
lines, which correspond to [b, a, X0, V mflX ] of [1.18", 2.84", -6.24", 3.37 km s _I ] and [1.43", 
4.57", 2.52", 5.27 km s -1 ], respectively. 












Table 1. 


Outflow parameters derived from HCN (3-2) lines 


Component 

Velocity interval 

(km s _1 ) 

V char 

(km S -1 ) 

R e// a 

(pc) 

t 

(10 3 yr) 

M 

(M 0 ) 

P 

(M 0 -km s~*) 

E 

(10 47 erg) 

M 

(10~ 4 M o yr~') 

L m 

(10 2 l q ) 

Low Blue 

[-30,-10] 

26.2 

0.1 

10.9 

12 

300 

0.8 

0.6 

0.6 

Low Red 

[30,40] 

29.5 

0.1 

7.8 

12 

350 

1.0 

0.9 

1.1 

High Blue 

[-80,-30] 

55.7 

0.2 

6.2 

6 

340 

1.9 

1.1 

2.6 

High Red 

[60,120] 

71.5 

0.2 

4.8 

10 

710 

5.0 

3.0 

8.7 


“corrected with an inclination angle of 20° with respect to the line of sight. R e ff is the full width at half-maximum (FWHM) 


deconvolved size. 
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